Optical Conductivity of Ferromagnetic Semiconductors 
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The dynamical mean field method is used to calculate the frequency and temperature depen- 
dent conductivity of dilute magnetic semiconductors. Characteristic qualitative features are found 
distinguishing weak, intermediate, and strong carrier-spin coupling and allowing quantitative deter- 
mination of important parameters defining the underlying ferromagnetic mechanism. 
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Observations of ferromagnetism in dilute magnetic 
semiconductors (DMS) 0-^] with transition tempera- 
tures as high as 110K (GaMnAs at 5% Mn (M]) or even 
above room temperature (GaMnN, GaMnP [H) have re- 
newed interest in these remarkable systems p- 1 1 ] , in part 
because of possible 'spintronic' applications |12|]. Several 
crucial scientific questions have emerged, most notably 
what parameters control the magnitude of the magnetic 
transition temperature and how these may be optimized. 
There is also an urgent need for experimental diagnos- 
tics, both for determining fundamental parameters such 
as the strength of the carrier-spin coupling, and for de- 
termining sample-specific parameters such as the carrier 
density and the degree of spin polarization of the mobile 
carriers. In this Letter we present theoretical calculations 
showing that measurements of the optical conductivity 
a(u),T) can be of great help in answering these ques- 
tions. Optical conductivity measurements have proven 
useful in understanding the physics of the colossal mag- 
netoresistance (CMR) manganites |D| where carrier-spin 
coupling is also crucial Q . 

We find that the DMS optical conductivity exhibits 
novel features not found in the CMR. To perform our cal- 
culations we use a recently developed non-perturbative 
method, the 'dynamical mean field theory' (DMFT) [D| 
to calculate a(u>,T) for the generally accepted model 
of dilute magnetic semiconductors. A non-perturbative 
method is needed because the crucial physics involves 
bound-state formation and other aspects of intermedi- 
ate carrier-spin couplings not accessible to perturbative 
methods. The system sizes required for direct numeri- 
cal simulations for this problem are impractically large, 
and obtaining accurate dynamical information from nu- 
meric is not easy. The DMFT method can handle dy- 
namics easily. Our results display striking, sometimes, 
counterintuitive, dependence of conductivity on carrier- 
spin coupling, carrier density, and temperature revealing 
key features of the underlying physics. Experimental ob- 
servation of our predictions should lead to crucial infor- 
mation about bound state formation and impurity band 
physics in this problem. 

It is generally believed [jyj that DMS are described by 
the generalized Kondo lattice model 

H = H host - ^ JSi ■ 4>t(Ri)&af34>f3(Ri) 
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where Hhost describes carrier propagation in the host 
semiconductor and the second (magnetic) term describes 
coupling of the carriers to an array of (impurity, e.g. Mn) 
spins at positions iij. The coupling has two sources: a 
spin-spin coupling and a potential scattering. Here we 
absorb the magnitude of the impurity spin into the cou- 
pling J (which we take to be positive), and represent 
the spin direction by the unit vector S. For simplicity 
we consider a host material with a single non-degenerate 
band; the extension to the multiband case relevant, for 
example, to the GaAs valance band is straightforward, in- 
volves no new features, and will be presented elsewhere. 
We therefore write 
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where Vr is a random potential arising 
magnetic defects in the material. 

The crucial physical issues are revealed by the con- 
sideration of a ferromagnetic state in which all impurity 
spins Si are aligned, say, in the z direction. Then the 
carriers with spin parallel to Si feel a potential — J + W 
on each magnetic impurity site and anti-parallel carriers 
feel a potential J + W. These potentials self-consistently 
rearrange the electronic structure. The spin-dependent 
part of this rearrangement provides the energy gain which 
stabilizes the ferromagnetic state. The key physics issue 
is, evidently, whether the potential W =F J is weak (so 
its effect on carriers near the lower band edge is simply 
a scattering phase shift) or strong (so only majority spin 
or perhaps both species of carriers are confined into spin- 
polarized impurity bands) . Recent density functional su- 
percell calculations [M suggest that in GaMnAs —J + W 
is close to the critical value for bound state formation for 
the majority spin systems. Then as temperature is in- 
creased, the Mn spins disorder and it is natural to ask 
how this physics changes. We shall see that all of this 
behavior is clearly revealed in the optical conductivity, 
which can therefore be used experimentally to probe the 
qualitative nature of the magnetic coupling. 

To compute the conductivity we write V — > (V — 
ieA/c) and apply the usual Kubo formula. To evalu- 
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ate the properties of H we employ the DMFT approx- 
imation, which amounts to assuming that the self en- 
ergy is momentum independent: S(p, lo) — > S(o;), so 
that the Green function G(p, oj) corresponding to H is 
G(p,w) = [u> — p 2 /2m — E(cj)] -1 . The self energy is 
given in terms of the solution of an impurity problem 
Jloj ]. The impurity problem requires a momentum cut- 
off, arising physically from the electron band-width. We 
impose the cutoff by assuming a semicircular density of 
states D(e) = \ftt 2 - e 2 /2irb with t = (2n) 2/3 /mb and b 3 
the volume per formula unit. The parameter t is chosen 
to correctly reproduce the band edge density of states. 
Within this approximation the real part of the conduc- 
tivity is given by 
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Our approximation for the density of states implies 
/ d 3 p /(2tt) 3 (p cos 9 /m) 2 -► / deD(e)<f>(e) with $(e) = 
(4* 2 -s 2 )/3@. 

The main panels of Fig. 1 show the evolution of 
the conductivity with temperature for two couplings; 
intermediate-weak (J = 1, Fig. 1(a)) and intermediate- 
strong (J = 1.5, Fig. 1(b)); the insets show the majority- 
spin densities of states. Consider first the T = T C curves 
(solid lines), where T c is the ferromagnetic transition 
temperature |11|. The J — 1 curve has approximately 
the Drude form expected for carriers scattering off ran- 
dom impurities (a closer examination reveals minor dif- 
ferences due to density of states variations near the band 
edge). In the J = 1.5 case the density of states plot shows 
that an impurity band is formed and the corresponding 
conductivity has two structures: a low- frequency quasi- 
Drude peak corresponding to motion within the impurity 
band and a higher frequency peak corresponding to ex- 
citations from the impurity band to the main band. We 
observe that the frequency of the upper peak does not 
directly give the separation between the impurity band 
Fermi level and the conduction band, because the vanish- 
ing of the velocity and density of states at the band edge 
means that the lowest unoccupied state is not optically 
active. 

We now turn to the temperature dependence. For the 
J = 1.5 curves we see that as T is decreased, the dc con- 
ductivity changes only very slightly, whereas the width 
of the low frequency peak increases and the high fre- 
quency peak moves up in energy and decreases in os- 
cillator strength. The increase in energy of the higher 
frequency peak may be understood from the density of 
states curves, which show a weak blue-shift of the conduc- 
tion band edge and a broadening of the impurity band. 
The counter-intuitive broadening (as T is decreased) of 
the lower frequency peak arises because, as the spins or- 
der, the binding of the carriers to the impurity spins 
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Temperature dependence of optical conduc- 
(a) J = 1.0*, x = 0.05, n = 0.02 with 



T c = 5.95 x 10-H/kB, and (b) J = 1.5*, x = 0.05, n = 0.02 
with T c = 1.24 x 10~ 2 t/kB- Insets show the density of states 
for both majority spin (thick curves) and minority spins (thin 
curves). The vertical lines show the zero temperature chemi- 
cal potentials. 

increases (as seen from the increase in separation between 
the chemical potential and the edge of the main band) 
corresponding to an increase in the basic scattering rate. 
The weak T-dependence of the dc conductivity occurs 
because the increase in scattering rate is compensated 
by another effect. Because the impurity band is spin- 
polarized carriers which are bound to impurity site must 
have spins parallel to impurity spin. Thus, as the spins 
order ferromagnetically, basic ability of carriers to move 
in the impurity band increases. This physics is famil- 
iar from the CMR materials [E3 and corresponds to an 
increase in conduction band oscillator strength. 

Consider next the J — 1 curve, where two effects occur 
as T decreases. First, the main quasi-Lorentzian peak 
decreases slightly in amplitude and increases slightly in 
half-width. The increase in width is due to increased 
carrier-spin coupling as mentioned above. Second, a new 
narrow peak appears. As can be seen from the inset in 
Fig. 1(a), at T — for this carrier concentration, the 
minority spin-band is slightly occupied and the sharply 
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FIG. 2. Imaginary part of the self energy as a function of 
the frequency for J = l.Ot and £ = 0.05. Solid (dashed) 
curves show the self energy at T = (T = T c ). 

peaked conductivity occurs because this small density of 
carriers is very weakly scattered. 

To understand the small scattering rate it is helpful 
to consider the Born approximation, which would lead 
to Im[E] ~ xN{Ef)J with x the magnetic impurity 
concentration, N(Ef) the density of states, and J the 
carrier-spin coupling. The combination of the small value 
of x and the small density of states at the band edge leads 
to a small scattering rate (in the Born approximation). 
Now consider corrections to the Born approximation. For 
the majority spin band, increase of J leads to the forma- 
tion of a spin-polarized bound state, so the corrections 
to the Born approximation must be such as to strongly 
increase the effective scattering rate. On the other hand, 
for minority spin carriers the increase of J leads to an an- 
tibound state at the top of the band, so that at the phys- 
ically relevant lower band edge, the corrections must be 
such as to decrease the effective scattering rate. Quanti- 
tatively, these effects are quite large. Fig. 2, for example, 
shows the calculated imaginary part of the T = major- 
ity and minority spin self energy for J = 1 . 

We now consider the density dependence of the conduc- 
tivity. Fig. 3(a) shows the evolution of a{uj,T — 0) with 
carrier density for J = 1. At very low density (n = 0.01) 
the minority spin band is unoccupied and the behavior 
associated with the majority spin band is observed. A 
sharp peak occurs when the minority band begins to be 
populated. Fig. 3(b) shows the evolution of <t(uj,T = 0) 
with n for the 'impurity-band' case J = 1.5. At very 
small n the conductivity is dominated by the 'impurity 
band' contribution, with a relatively weak feature corre- 
sponding to transitions from the impurity band to the 
main band. As the carrier density is increased the band 
oscillator strength in transitions between the impurity 
and main band increases dramatically, both in absolute 
terms and relative to the intra- impurity band transitions. 
When the Fermi level crosses into the minority spin band 
(note that the minority spin impurity band is at the up- 
per band edge; here we have only main-band states) an 
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FIG. 3. Density dependence of optical conductivity for (a) 
J = l.Ot, x = 0.05, and (b) J = 1.5t, x = 0.05 for T = OK. 
Inset shows the density of states and the chemical potentials 
at zero temperature. For n > n c the minority subband is 
occupied. 

additional sharp low frequency component is observed. 

Fig. 4 shows the sensitivity of the predicted behavior 
to potential scattering. The main panel shows the evo- 
lution of the conductivity as the scalar part W of the 
electron-impurity potential is varied. We study here a 
case in which at zero scalar potential the impurity band 
is well formed. As the potential is made more attractive, 
the impurity band features become more pronounced; as 
it is made more repulsive, the impurity band rapidly re- 
joins the main band and the extra feature in the con- 
ductivity is lost. The inset shows the dependence of the 
transition temperature on density at different potential 
strengths. We see that for attractive potential the depen- 
dence on potential strength is weak, but as the potential 
is made more repulsive the impurity band is destroyed 
and T c decreases. 

It is instructive to compare the optical proper- 
ties of DMS to those of another system with strong 
carrier-spin couplings, namely CMR manganites such as 
(Lai_ :E Ca x )Mn03 [@. In CMR, instead of being dilute 
random impurities as in DMS, the Mn ions form an or- 
dered lattice. They possess large local moment, to which 
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FIG. 4. Optical conductivities for various of impurity po- 
tential scattering. Inset shows the transition temperature as 
a function of density. 

mobile carriers are very strongly coupled. Thus instead of 
a spin-polarized impurity band, there is a spin-polarized 
conduction band, sufficiently well separated from other 
bands in the solid that the contribution to the optical 
conductivity arising from it may be clearly identified and 
analyzed |L3|] . The periodic arrangement of the Mn sites 
means that (in the absence of other physics) the scatter- 
ing rate decreases as T is lowered unlike in the DMS sys- 
tem. The conductivity turns out to be most usefully char- 
acterized by its integrated area (spectral weight), and the 
change of this quantity was shown to be a good predictor 
of the magnetic transition temperature 

These CMR ideas have limited applicability to the 
DMS systems. For J < 1 the physics is of extended states 
scattered by defects, and this is sufficiently different from 
the CMR situation that the spectral weight in DMS is not 
related in a useful manner to T c in this weak coupling sit- 
uation. However, as J is increased and the spin-polarized 
impurity band develops the physics becomes more anal- 
ogous to that of CMR. For example, in Fig. 1(b) the 
two structures in the conductivity are relatively well sep- 
arated, and the spectral weight in the lower feature in- 
creases as T decreases below T c . For J > 2 (not shown) 
the low feature becomes completely separated from the 
upper one. We therefore define an effective 'impurity 
band spectral weight' for J > 1 by integrating the con- 
ductivity from to = to the conductivity minimum. We 
find again that the changes in 'impurity band spectral 
weight' are very weak in DMS systems relative to those 
observed in CMR materials, so that spectral weight ideas 
useful in CMR do not carry over to DMS. This is an im- 
portant qualitative distinction between DMS and CMR 
materials. 

To summarize, we have shown that the frequency, den- 
sity, and temperature dependence of the conductivity 
contains important information about the physics of di- 
lute magnetic semiconductors; in particular the forma- 
tion of a spin-polarized impurity band leads to a peak 



centered at a non-zero frequency. We present several, at 
first sight, counterintuitive findings. We find an increase 
in scattering rate as T is decreased, signalling enhanced 
carrier-spin coupling with increasing spin alignment. We 
also show that in certain doping and coupling regimes a 
very narrow conductivity peak could occur correspond- 
ing to a slightly occupied minority spin band. It occurs 
because the repulsive interaction between local moments 
and "wrong-spin" carriers suppresses the carrier ampli- 
tude at the impurity site, reducing the effective carrier- 
spin coupling and consequently narrowing the wrong-spin 
conductivity peak; in addition in three dimensions the 
vanishing of the density of states at the band edge fur- 
ther decreases the rate. This suggests that a small oc- 
cupation of the minority band could be quite dangerous 
from the spintronic applications point of view, particu- 
larly in three dimensional devices. 

Finally, we note that the most interesting phenom- 
ena involve intermediate couplings, intermediate temper- 
atures and non-zero frequency response. This regime is 
very difficult to treat by standard analytical or numeri- 
cal methods; it is therefore fortunate that the dynamical 
mean field method allows access to this regime. 
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